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Abstract. A nucleation picture of magnetization switching in single-domain 
ferromagnetic nanoparticles with high local anisotropy is discussed. Rele- 
vant aspects of nucleation theory are presented, stressing the effects of the 
particle size on the switching dynamics. The theory is illustrated by Monte 
Carlo simulations and compared with experiments on single particles. 



1. Introduction 

The dynamics of magnetization switching in nanometer-sized particles of 
highly anisotropic ferromagnets is interesting, from both the scientific and 
technological points of view. The basic scientist sees in such particles a lab- 
oratory to study the decay of a metastable phase towards equilibrium, while 
the technologist sees a promising material for ultrahigh-density magnetic 
recording media. Although ferromagnetic nanoparticles have been studied 
experimentally for a long time [1], until recently this was only possible 
with powders. However, with modern techniques of nanofabrication [2] and 
ultrahigh-resolution methods to detect the magnetization, such as Magnetic 
Force Microscopy (MFM) [3] Lorentz microscopy [4] and micro-SQUID de- 
vices [5] , one can now synthesize and study such particles individually. 
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The most common description of magnetization switching is a mean- 
field approach, originally due to Neel [6] and Brown [7]. To avoid an energy 
barrier due to exchange interactions of strength J, uniform rotation of all 
the atomic moments in the particle is assumed. The remaining energy bar- 
rier, A, is caused by magnetic anisotropy, which is a combination of crystal- 
field and magnetostatic effects. The equilibrium thickness of a wall between 
oppositely magnetized domains is £ oc \J J/ A. For particles smaller than £, 
the uniform-rotation picture is reasonable. If the anisotropy is largely mag- 
netostatic, the resulting demagnetizing field causes particles larger than £ 
to form oppositely magnetized domains, and switching is achieved through 
the field-driven motion of preexisting domain walls. However, if the an- 
isotropy is largely due to the local crystalline environment, there exists a 
window of particle sizes that are larger than £ but smaller than the size 
at which the particle becomes multidomain. [This is for instance often the 
case in ultrathin films.] Such particles can be modeled as Ising systems with 
local spin variables, S{ = ±1. Depending on the degree of anisotropy, these 
spins can either represent the z component of individual atomic moments, 
or one can coarse grain the system by rescaling all lengths in terms of £, so 
that the represent block spins. The Ising Hamiltonian is 

H = -J^SiSj - H^Si . (1) 

(i,J> * 

Here J > is the ferromagnetic exchange interaction, H is the applied 
magnetic field times the local magnetic moment, and the sums J2{i,j) an d 
J2i r un over all nearest-neighbor pairs and all sites on a suitable lattice, re- 
spectively. Here we only report numerical results for two-dimensional square 
lattices, but our theoretical arguments are valid for general spatial dimen- 
sion. The order parameter is the dimensionless magnetization, 

m = N- 1 Y,s l , (2) 

i 

where N is the total number of Ising spins in the particle. 

In the highly anisotropic nanoscale ferromagnets described by Eq. (1) 
[and modifications discussed below], the state of uniform magnetization op- 
posite to the applied field is properly viewed as a metastable phase. This 
nonequilibrium phase decays by neither uniform rotation nor by the motion 
of preexisting domain walls, but rather by the thermal nucleation and sub- 
sequent growth of localized droplets, inside which the magnetization is par- 
allel with the field [8]. This decay mechanism yields results very similar to 
effects observed in recent experiments on well-characterized single-domain 
ferromagnets in the nanometer range. In this paper we concentrate on a 
maximum in the switching field (or coercivity) vs. particle size [3]. Another 
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quantity which is often measured in experiments, is the probability that 
the particle has not switched within a specified waiting time [4, 5]. Results 
concerning this quantity can be found in Refs. [8] and [9]. 

The Ising model does not have an intrinsic dynamic. To simulate the 
effects of thermal fluctuations we therefore use a local stochastic dynamic 
which does not conserve the order parameter, such as the ones proposed by 
Metropolis et al. [10] or Glauber [11]. In order to perform simulations on 
the very long timescales necessary to observe metastable decay, we use a 
so-called "rejection- free" Monte Carlo (MC) algorithm [12]. The basic time 
scale of the MC simulation [MC Steps per Spin (MCSS)] is not known from 
first principles and must be fitted to experiments. It is expected to be on 
the order of a typical inverse phonon frequency, 10 _9 -10 -13 s. 



2. Nucleation and Growth 

This section is a brief primer on the theory of nucleation and growth as it 
applies to systems in the dynamic universality class of kinetic Ising models 
with nonconserved order parameter. For further details, see Refs. [8, 13, 14, 
15, 16]. 



2.1. BACKGROUND 

The central problems in nucleation theory are to identify the fluctuations 
that lead to the decay of the metastable phase and to obtain their free- 
energy cost, relative to the metastable phase. For Ising- like systems with 
short-range interactions, these fluctuations are compact droplets of radius 
R. The magnetization inside the droplet is parallel with the applied field and 
has a magnitude near the temperature dependent zero-field magnetization, 
m sp (T), which is nonzero below the critical temperature, T c . The free energy 
of the droplet has two competing terms: a positive surface term oc R d ~ l , 
and a negative bulk term oc \H\R d , where d is the spatial dimension. The 
competition between these terms yields a critical droplet radius, 

(d-l)a(T) 

where a(T) is the surface tension. Droplets with R < R c most likely decay, 
whereas droplets with R > R c most likely grow further to complete the 
switching process. The free-energy cost of the critical droplet (R = R c ) is 

A FsM HX> = aMTf(^^) d \ (4) 

where is a weakly T dependent shape factor such that the volume of 
a droplet of radius R equals Q.<iR d . The subscript SD stands for Single 
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Droplet, as explained below. Nucleation is a stochastic process, and the 
nucleation rate per unit volume is given by a Van't Hoff-Arrhenius relation: 

I{H,T) oc |#|^exp 

where k B is Boltzmann's constant, and S(T) is the iif-independent part of 
AFsd- The prefactor exponent K equals 3 for the two-dimensional Ising 
model and —1/3 for the three-dimensional Ising model [16, 17, 18]. 



AF SD (H,T) 
k B T . 



\H\ K exp 



k B T\H\ d - 



, (5) 



2.2. EFFECTS OF FINITE PARTICLE SIZE 

For particles of finite linear size, L, an important crossover occurs for combi- 
nations of H, T, and L, such that R c « L. This yields a T and L dependent 
crossover field called the Thermodynamic Spinodal (ThSp) [15, 16]: 

(d-lMT) 

i7 ThSp (T,L)^ . (6) 

For < i^ThSp) -Rc would exceed L. This is called the Coexistence 
(CE) regime because the critical fluctuation in such weak fields resem- 
bles two coexisting slabs of opposite magnetization [15, 16]. The average 
metastable lifetime in the CE regime is approximately 



tce(H,T,L) ~ exp 



2a(T)L d - 1 - 2Am sp {T)\H\L d 
kB~T 



(7) 



where A is a nonuniversal constant. Since \H\ < -f/rhSp ~ L^ 1 , the domi- 
nant size dependence is an exponential increase with L d—1 . This behavior 
also holds for more general boundary conditions than the periodic boundary 
conditions used to obtain Eq. (7) [14]. 

For \H\ > ifrhSp (but not too large, as we shall see below), the lifetime 
is dominated by the inverse of the total nucleation rate, 



t sd (#, T, L) w (L d I(H, T)) 1 oc L- d \H\ K exp 



S(T) 



k B T\H\ 



d-i 



(8) 



It is inversely proportional to the particle volume, L d . The subscript SD 
stands for Single Droplet and indicates that in this regime the switching is 
completed by the first droplet whose radius exceeds R c . 

A second crossover, called the Dynamic Spinodal (DSp) [15, 16], is pre- 
dicted when one observes that a supercritical droplet grows at a finite 
velocity, which for large droplets is proportional to the field: v ~ v\H\. A 
reasonable criterion to locate the DSp is that the average time between 
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nucleation events, tsd, should equal the time it takes a droplet to grow to 
a size comparable to L. This leads to the asymptotic relation 



H DSp (T,L) 



(d-1) 
2m sp (T) 



Q d a(TY 



(d+l)k B TlnL 



i 

d-1 



(9) 



For \H\ > i^DSp) the metastable phase decays through many droplets 
which nucleate and grow independently in different parts of the system. 
This is called the Multidroplet (MD) regime [15, 16]. A classical theory 
of metastable decay in large systems [19, 20, 21] gives the lifetime in this 
regime, 

' I{H,T)n d {u\H\f 



t M d(H,T) 



(d + l)ln2 



d+l 



(10) 



which is independent of L. 

The switching field, H sw (t w ,T, L), is the field required to observe a 
specified average lifetime, t w . It is found by solving Eqs. (7), (8), and (10) 
for H with t w for the respective average lifetimes, tce, etc. The resulting 
L dependence of H sw is illustrated by the MC data shown in Fig. 1(a). 
It consists of a steep increase with L in the CE regime, peaking at the 
ThSp, followed by a decrease in the SD regime towards a plateau in the 
MD regime. 



3. Numerical Results 

In this section we present some representative results of simulations of 
two-dimensional Ising ferromagnets, which we compare with the theoretical 
predictions of the previous section and with experiments. 



3.1. PURE SYSTEM WITH PERIODIC BOUNDARY CONDITIONS 

The simplest model considered is a two-dimensional square-lattice Ising 
system with periodic boundary conditions. The switching fields for this 
model at T = 0.8 T c are shown in Fig. 1(a) for t w = 100 and 1000 MCSS. 
The L and t w dependencies expected from the results of Sec. 2.2 are clearly 
seen. We emphasize that the decrease in the SD region is not due to an 
equilibrium domain structure. It is an entropy effect of purely dynamical 
origin, arising from the volume factor in Eq. (8) [14] . Analogous corrections 
to nucleation rates in fluids were proposed by Lothe and Pound [22]. 

For qualitative comparison we show in Fig. 1(b) effective switching fields 
for nanoscale Ba-ferrite particles, obtained by MFM experiments [3]. We 
propose that the peak observed in the switching field may be of the same 
purely dynamical origin as in kinetic Ising models. 
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Figure 1. Switching fields vs. particle size, (a): MC simulations for a two-dimensional 
Ising ferromagnet with periodic boundary conditions. The dotted line is the ThSp, and 
the dashed line is the DSp. After Rcf. [8]. (b): Effective switching fields for nanoscale 
Ba-ferrite particles. Data digitized from Fig. 5 of Ref. [3] . 

3.2. EFFECTS OF A DEMAGNETIZING FIELD 

A reasonable objection to the model defined by Eq. (1) is the absence of 
dipolar interactions, which causes it to be single domain for all L. To address 
this shortcoming without the large computational expense of recalculating 
dipole sums at every step in the dynamical simulation, a model was in- 
troduced in which the demagnetizing field was approximated by adding a 
weak long-range antiferromagnetic term: Hd = Ho + DL d m 2 [13]. Particles 
smaller than Lo ~ 2a(T)/D remain single domain [23], but the demagne- 
tizing factor D decreases the free-energy barrier towards nucleation of the 
equilibrium phase. Addition of the demagnetizing factor was found to re- 
duce the average lifetime by an analytically predictable amount, as shown 
in Fig. 2. However, no qualitative differences from the behavior described 
above were observed. 

3.3. HETEROGENEOUS NUCLEATION 

Next we discuss ways in which the homogeneous nucleation observed in pure 
systems with periodic boundary conditions is modified by heterogeneous 
nucleation at the particle surface or at quenched inhomogeneities. 

3.3.1. Modified Boundary Conditions 

The use of periodic boundary conditions allows one to study bulk nucleation 
without complications due to the particle surface. Since the surface can be 
modified in various ways by reconstruction, adsorption, oxidation, etc., one 
cannot in general predict whether it will enhance or inhibit nucleation. 
However, even addition to the Ising model of a surface field or modified ex- 
change interactions at the surface produces complicated crossovers between 
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Figure 2. Relative changes in the average metastable lifetime versus the reduced demag- 
netizing factor, x — 2Dm sp (T)/\H\, at T = 0.8T C . The solid curves are analytical results 
that only require parameters determined for D = 0. After Ref. [13]. (a): SD regime, 
\H\ = 0.2J, L = 10. (b): MD regime, L = 100. 



surface and bulk nucleation [14]. In general, the changes reduce the height 
of the peak in H sw vs. L, but for a wide range of modifications it remains 
clearly discernible. Examples are shown in Fig. 3(a). 

3.3.2. Quenched Randomness 

Another way in which heterogeneous nucleation may dominate, is through 
quenched impurities. An exploratory study was presented in Ref. [24]. Bond 
dilution was observed to reduce H sw by a factor approximately independent 
of L, as shown in Fig. 3(b), while random spin magnitudes led to non-self- 
averaging behavior and a wide distribution of lifetimes. 

3.3.3. Coercivity of Fe Sesquilayers on W(110) 

Much interest has recently been devoted to ultrathin iron films on W(110) 
substrates [25, 26, 27, 28, 29]. The so-called sesquilayer systems, which 
consist of islands of a second monolayer of Fe on top of an almost perfect 
first monolayer [27], have particularly interesting magnetic properties [28, 
29, 30]. Around an Fe coverage of approximately 1.5 monolayers (ML), the 
coercivity exceeds that of a monolayer or a bilayer by more than an order 
of magnitude [29]. 

Magnetization switching in this system is expected to occur through 
the field-driven motion of preexisting domain walls, which are pinned at 
the second-layer islands. Based on this picture, the coercivity has been 
calculated by micromagnetic methods [28, 29]. However, those calculations 
did not consider thermal effects and were also essentially static. 

To account for thermal depinning and the dependence of the coercivity 
on the frequency of the applied field, a two-layer Ising model has been de- 
veloped for this system [31]. This is a reasonable approximation since the 
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Figure 3. Effects of heterogeneous nucleation on H sw at T ~ 0.57 T c with 
t w = 30000 MCSS. For comparison, the top curve in both panels corresponds to ho- 
mogeneous nucleation in a pure system with periodic boundary conditions, (a): Effects 
of boundary conditions in a pure system. Middle curve: square system with periodic 
boundary conditions in one direction and free boundary conditions in the other. Bottom 
curve: circular system with free boundary conditions. Data from Ref. [14] . (b) : Effects of 
random bond dilution in a system with periodic boundary conditions. After Ref. [24] . 



crystal-field anisotropy for Fe monolayers on W(110) is almost two orders 
of magnitude larger than for bulk Fe [25]. Simulations were performed on 
a computational lattice in which the second-layer island morphology was 
reproduced using STM images of real systems [27], and the exchange in- 
teractions were chosen to reproduce the experimentally observed critical 
temperature of an Fe monolayer, T c = 230 K [26]. 

Two simulation snapshots of the magnetic domain wall moving across a 
sesquilayer are shown in Figs. 4(a) and (b). The wall moves intermittently, 
spending most of its time pinned against the "windward" side of the islands. 
Thermally nucleated depinning events are followed by rapid advances to 
the next metastable position. The coercivity is estimated from the average 
domain-wall velocity. Estimated coercivities for two different temperatures 
and driving frequencies are shown vs. the Fe coverage in Fig. 4(c). The 
experimentally observed nonmonotonic coverage dependence [29] is repro- 
duced, as well as the temperature [28] and frequency [30] dependencies. The 
model yields an approximately linear dependence of the inverse coercivity 
on the logarithm of the frequency, shown in Fig. 4(d). Over a few decades 
of frequency, this is hard to distinguish numerically from a power law [30] . 

4. Conclusions 

We have presented a brief overview of a nucleation theory of magnetiza- 
tion switching in single-domain ferromagnets in the nanometer range. We 
emphasized the dependence of the switching field or coercivity on the par- 
ticle size and demonstrated that the model is capable of reproducing the 
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Figure 4- Simulations of magnetization switching in Fe sesquilayers on W(110). After 
Ref. [31]. (a) and (b): Snapshots of a domain wall propagating across a sesquilayer with 
coverage 1.26 ML. The high-contrast region represents the growing equilibrium phase. 
The area shown is 654 A x 610 A [109 x 102 computational cells], and the island con- 
figuration was digitized from Fig. 1 j) of Ref. [27]. The simulated temperature and field 
correspond to 132 K and 0.26 T, respectively. The time elapsed between the two snapshots 
is approximately 1.5xl0 6 MCSS, corresponding to 1.5xl0~ 6 s. A movie of this simula- 
tion is found at http://www.scri.fsu.edu/~rikvold. (c): Sesquilayer coercivity vs. Fe 
coverage, estimated by extrapolation to weak fields. The lower curve should be compared 
with Fig. 3 of Ref. [29]. (d): The frequency dependence of the estimated coercivity. 

experimentally observed maximum in the switching field vs. particle size. 

Our discussion places the switching dynamics of nanoscale ferromagnets 
in the context of metastable decay in finite systems. This interdisciplinary 
field is experiencing a renaissance due to new methods of nanofabrication 
and observation of individual systems. In addition to magnets, results have 
recently been published for systems as different as liquid mixtures [32] and 
semiconductor nanocrystals [33]. 
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